Fermi super- Tonks-Girardeau state for attractive Fermi gases in an optical lattice 
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We investigate properties of an excited state of strongly attractive Fermi Hubbard model, in 
analogue with the super- Tonks-Girardeau state. In contrast to the ground state of the attractive 
Hubbard model, such a state is the lowest scattering state with no pairing between attractive 
fermions and can be viewed as a Fermi generalization of super- Tonks-Girardeau state in the optical 
lattice system. We show that this state can be realized in the spin-1/2 Fermi optical lattice system by 
a sudden switch of interaction from the strongly repulsive regime to the strongly attractive regime. 
With the aid of Bethe-ansatz method, we calculate energies of both the Fermi Tonks-Girardeau 
gas and the Fermi super- Tonks-Girardeau state of spin-1/2 ultracold fermions and show that both 
energies approach to the same limit as the strength of the interaction goes to infinity. By exactly 
solving the quench dynamics of the Hubbard model, we demonstrate that the Fermi super- Tonks- 
Girardeau state can be realized efficiently. This allows the experimental study of properties of Fermi 
super- Tonks-Girardeau gas in optical lattices. 

PACS numbers: 37.10.Jk, 03.75.Ss, 03.75.Kk 



I. INTRODUCTION 

Ultracold atomic gases trapped in one-dimensional 
(ID) waveguides have become one of the most active re- 
search field of cold atom physics in recent years [l| . Due 
to their good tunability, the ID atomic gases have pro- 
vided an ideal platform for studying and testing the the- 
ory of low-dimensional many-body systems 0,|3[. Tuning 
the effective interaction strength between atoms via Fes- 
hbach resonance has led to the experimental realization 
of Tonks-Girardeau (TG) gases [Hy], which describes the 
Bose gas in the strongly repulsive limit and exhibits the 
feature of fermionization. By switching the interaction 
between atoms of an initial TG gas from strongly repul- 
sive to strongly attractive, the experimental observation 
of a ID super Tonks-Girardeau (STG) gas [3] of bosonic 
cesium atoms was reported very recently [f| ■ In contrast 
with the strongly repulsive interacting TG gas, bosons in 
STG state interact via strongly attractive interaction. A 
surprising feature of this many-body state is its good sta- 
bility even under strongly attractive interaction instead 
of decaying into the lower atomic bound states . The 
stability of the STG gas could be well understood from 
the quench dynamics of the ID integrable Bose gas @. 

The experimental realization of stable excited quan- 
tum gas phase opens a new area for searching exotic 
quantum phases in ultra-cold systems [9l4l9|, Particu- 
larly, the exotic experimental results |5| have stimulated 
intensive theoretical studies of the STG gases from vari- 
ous aspects So far, most of the theoretical works 
on the STG gases have focused on the bosonic gases in 
continuum systems. In this work, we study the possible 
realization of the Fermi super- Tonks-Girardeau (FSTG) 



'Electronic address: schen@aphy.iphy.ac.cn 



gas for a Fermi gas loaded into a deep ID optical lattice, 
which is described by the Fermi Hubbard model [2(| • The 
FSTG state for a ID Fermi continuum gas was studied 
in Ref. [IH, The current study can be viewed as 

a generalization for the realization of FSTG gas in the 
lattice systems. Despite that in the low-density limit the 
Hubbard model shares similar behaviors with the con- 
tinuum Yang-Gaudin model [2ll[23 |. the Hubbard model 
exhibits some new features due to the existence of the 
band structure for the lattice model. One of the new 
features is the existence of repulsive bound pairs in the 
high bands which is absent in the continuum model. The 
other one is the existence of the Mott insulating phase in 
the half-filling case. Although the Hubbard model is one 
of the fundamental model in condensed matter physics, 
most of previous theoretical works focused on its ground 
state and thermodynamical properties. Our study shall 
shed lights on properties of some highly excited states 
which can be accessible in current experimental condi- 
tions. 

Stimulated by the experiment of the bosonic STG gas, 
we first prepare an ultracold Fermi gas in the strongly 
repulsive regime, and then suddenly switch the inter- 
action to the strongly attractive regime. Through this 
way, we can reach a stable highly excited phase which 
is the lowest scattering state of the attractive Fermi gas. 
To understand properties of the STG state of the at- 
tractive Hubbard model, we shall analyze the spectrum 
structure of the Hubbard model for both the repulsive 
and attractive cases. By calculating the energy of the 
FSTG state analytically based on the Bethe-ansatz (BA) 
method, we show that the energy of FSTG gas state in 
the strongly attractive limit approaches the same limit of 
the ground state energy of ID strongly repulsive Fermi 
Hubbard model. This implies that the FSTG state can 
be accessible from the the ground state of the strongly 
repulsive Fermi gas by a sudden switch of interactions, 
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which is also confirmed by the exact calculation of the 
quench dynamics of the ID Fermi gas on the optical lat- 
tice by using numerical exact diagonalization method. 



II. MODEL AND FSTG STATE 

We consider a ID ultracold Fermi gas composed of 
N = Nf + iVj. spin-1/2 fermionic atoms in a deep optical 
lattice, which can be well described by the well-known 
Hubbard model (HM), 



H=-t 



i,<r 



U J2 



n it n. 4 , (1) 



where c\ a {ci^) with a =f, \, is the creation (annihilation) 
operator of fermions at the ith site, t and U denote the 
hopping amplitude and the on-site interaction strength, 
respectively. The ratio U/t can be tuned by varying the 
depth of the optical lattice potential and by using Fesh- 
bach resonance technique. For convenience, we set t = 1 
as the energy scale. Without loss of generalization, we 
assume that N\, < Nf. 

The one-dimensional Hubbard model ([T]) with peri- 
odic boundary condition is exactly solvable by the Bethe- 
ansztz (BA) method [2(| with the BA wavefuction 

(f>(xi, ...,x N ) =^2^2Q(xqi < ■■■ < xqn)x 
Q p 

N 

[Q,P]exp[i^2k P jX Qj ], (2) 

3=1 

where kiS represent quasimomenta, Ps and Qs repre- 
sent permutations of fc^s and ijs, respectively. For the 
eigenstate with the total spin S = N/2 ~ M (M = 
Ni), the coefficient [Q,P] can be explicitly expressed as 

[Q,P] = e(Q 1 )e(Q 2 )6(y 1 ,y 2 ,...,y M ) i &(yi,y2,-,yM;P), 
where Q\ is the ordering of the first M fermions 
and Q2 the ordering of the rest fermions, y±, y%, yu 
are the coordinates of down spins in the lat- 
tice with length L, and ^{yx, yi-, yM'i P) = 

z R <r) n K * ( A ^- - A «< - w / 2t ) n,=i wxi 1 - 

A Rl + iU/At) Ut= yi +i ( sin k Pt - Affi - iU/At)}. The pa- 
rameters kjS and A a s are determined by the Bethe-ansatz 
equations (BAEs) [20| : 
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where L is the size of the optical lattice. The eigenvalues 



are given by E = — 2t^2j—% cos kj. The structure of the 
solution of BAEs of Hubbard model is relevant to the 
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FIG. 1: (Color online) Comparison between exact solutions 
(denoted by O) an< i approximate solutions (denoted by +) 
of BAEs. (a) Quasi-momentum distributions for the ground 
state of the repulsive Fermi gas and the FSTG state of the 
attractive Fermi gas with different values of U/t and t is set 
to 1 as energy scale here, (b) The corresponding solutions of 
A Q for different values of U/t. 



filling factor of N/L. In the following, we shall consider 
the case with N/L < I. 

The BAEs (J3J and (gj hold true for both the repulsive 
and attractive U, however the structure of the solutions 
is quite different for U > or U < 0. For the repul- 
sive interaction with U > 0, both the solutions of kj 
and A a for the ground state (GS) and low excited states 
are real numbers. The ground state solution corresponds 
to Ij = (N + l)/2 - j and J a = (M + l)/2 - a. In the 
strongly repulsive interaction limit U/t — > 00, the ground 
state energy is identical to that of a polarized N-fermion 
system j23[. On the other hand, if the on-site interac- 
tion between fermions with different spin is attractive, 
i.e. U < 0, the ground state is then composed of N — 2M 
real kiS and 2M complex ones. In the strongly attrac- 
tive interaction limit U/t <C 1, the complex solutions 
take the 2-string form [24]]: sink a ss A a + i\U\/At, and 
s'mkM+a ~ Aq, — i\U\/At. Besides the complex solutions, 
the BAEs also have real solutions for U < 0, which de- 
scribe the scattering states of attractive fermions. The 
Fermi super- Tonks- Girardeau (FSTG) gas state corre- 
sponds to the lowest real solutions of BAEs ([3]) and ^ 
with U < 0. These scattering states are gas-like excited 
states of the attractive spin-1/2 ultracold fermions which 
are above states including at least one paired bound 
state, while the ground state of the system is composed 
of M tightly bound fermion pairs. 

Next we explore the scattering solution of the Bethe- 
ansatz equations in the strongly interacting limit. As 
\U\/t 00, the solution of A Q is proportional to U, 
whereas sin kj is always finite with | sin kA < 1 . There- 
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fore the quasimomenta can be given approximately 



kjL = 2-kIj+Aq 



where 
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with U 1 — U /At. Here we consider the case with Nf = 
N±. Under this condition, the values of A Q are symmetric 



about zero. It follows that Aq 
ingly Eq.Q is simplified as 



Ai = and correspond- 
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where <r = A\. In general, Ax(U) ^ Ax(—U) since 
the solution A Q of Eq. (0| are not symmetric for U 
and —U. However, in the strongly interacting limit, up 
to the order of J7" 1 Eq. (gj) becomes 2 tan" 1 (A^) = 

jf2nJ a + jj J^pLx 2 tan" 1 ( A J j which has the same 
form as BAE of the Heisenberg spin chain and is invariant 
under the operation P : U — > — U, A a — > — A a . Therefore, 
we have Ai(?7) = Ax(-U) up to the order of C/" 2 . It 
follows that the ground state energy of the Fermi Tonks- 
Girardeau (FTG) gas in the strongly repulsive limit and 
the energy of FSTG state in the strongly attractive limit 
are given by 
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where I, = (N + l)/2 - j for both the FTG and FSTG 
gas. Here, for convenience, we call the ground state of 
the spin- 1/2 Fermi gas in the strongly repulsive limit as 
the FTG state. Obviously, in the limit of \U\ — > oo, we 
have Efstg = Eftg- In Fig- HI we make a comparison 
between the exact solutions given by numerically solving 
the Eq. ([3|) and (gj directly and the approximate solu- 
tions given by solving Eq. © iteratively. For an example 
systems with L = 100,A^=10 and M = 5, we show that 
for large enough values of U/t, the approximate solutions 
agrees very well with the exact solutions. In Fig. [T](a), 
one can see that the quasimomentum distributions for the 
ground state of repulsive Hubbard model and the FSTG 
state approach the same limit from different sides when 
\U\/t goes infinite. Correspondingly, Efstg and Eftg 
also approach the same limit as shown in Fig. [51 




|U| 



FIG. 2: (Color online) The energies Eftg (solid line) and 
Efstg (dashed line) vs U. Dots in the figure denote the 
approximate solutions via expansion in strongly interacting 
limit. For large enough \U\, they agree very well with the 
exact solutions. 



III. PREPARATION OF FSTG STATE 

The FSTG state can be realized in a ID deep optical 
lattice by a sudden switch of interaction similar to the 
experimental realization of bosonic STG gas in Ref. Q ■ 
Suppose that the initial state \^ini(t = 0)) = \ipo(Uo)) 
is prepared at the ground state in the strongly repulsive 
regime with Uo/t 3> 1, after a sudden switch to the op- 
posite regime with interaction strength U/t <C —1, the 
wave-function = e ~ iH( - u ^ \^ in i{U )) can be calcu- 

lated via 



(7) 



where c„ — (ip n (U) \ ipo(Uo)) with \ip n (U)) representing 
the n-th eigenstatc of the Hubbard model with on-site 
interaction U . It is straightforward that |c„| 2 is the tran- 
sition probability from the initial state to the n-th eigen- 
state of H(U). 

To give a concrete example which may help us get an 
intuitive understanding of the properties of the FSTG 
state of the attractive Hubbard model, we display the 
full energy-momentum spectra of the Hubbard model 
with L = 30, N = 4, M = 2, and U = ±15 in Fig. 
[3] As shown, the spectra is spilt into a series of sep- 
arated bands. For the repulsive case, the lowest band 
is a scattering continuum of N (here N = 4) unpaired 
fermions. The middle band is a scattering continuum 
formed by one tightly bound fermion pair and two un- 
paired fermions, whereas the top band is a scattering 
continuum of two tightly bound fermion pairs. The spec- 
tra for the attractive case is similar but in reverse order. 
The gap between centers of neighboring bands equals ap- 
proximately to the binding energy \U\ of a fermion pair. 
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FIG. 3: Full spectra of energies vs. total momentum K 
for the HM with L = 30, N = 4, M = 2, U = 15 (left) and 
U — —15 (right). The hopping term t is set to 1. 



FIG. 5: Doubly occupied sites Nd versus the interaction U 
for a system with L = 20, N = 10, M = 5. (a) STG state for 
system with U < and ground state for system with U > 0. 
(b) Ground state for system with U < 0. 
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FIG. 4: (a) The transition probability from the initial ground 
state of the repulsive system with Uo = 100 to the FSTG state 
after the sudden interaction switch to the attractive regime, 
(b) The transition probability from the initial ground state of 
the attractive system with Uo = —100 to the highest excited 
state of the top band of repulsive Hubbard model after the 
sudden switch of interaction. 



These separated bands are no longer distinguishable as 
the interaction strength is comparable to the bandwidth. 
The zero-momentum lowest scattering state of TV un- 
paired fermions denoted by a red star in Fig. EJa) is 
just the ground state of the repulsive Fermi gas, whereas 
the red star in Fig. [3jb) indicates the Fermi super- Tonks- 
Girardeau gas state. 

Starting from the ground state of the HM with a re- 
pulsive interaction and then suddenly switching the in- 
teraction to the attractive side, we evaluate the transi- 
tion probabilities from the initial repulsively GS to each 



eigenstate of the attractive Hubbard model. As shown 
in Fig. [4^,, we find that the transition probability to the 
lowest state (FSTG) of a given S in the top scattering 
band is very close to 1, when both the repulsive interac- 
tion of the initial state and the attractive interaction of 
the final state are strong enough. As the transition prob- 
ability to the lowest scattering phase is almost 1 in the 
strongly interacting regime, the transition probability to 
the lower paired states is almost completely suppressed, 
thus we expect that such a highly excited gas-like state 
of the strong attractive Fermi gas in optical lattice can 
be experimentally realized. Actually, the stable excited 
scattering state prepared in this way can be viewed as a 
Fermi generalization of the STG gas in the optical lattice. 
If the system enters to the weakly interacting regime, 
the transition probability to the FSTG state decreases 
quickly, whereas the transition probability to the ground 
state increases. 

Due to the existence of band structure, there exist 
states of repulsively bound pairs above the lowest contin- 
uum band for the repulsive Hubbard model. Such kind of 
repulsively pairing state is absent in the continuum Yang- 
Gaudin system [131 ]. Particularly, the top band of the 
strongly attractive Hubbard model is completely com- 
posed of repulsively bound pairs, which is very similar to 
the ground state of attractive Hubbard model composed 
of attractively bound pairs. Next we show that the high- 
est excited state composed of repulsive bound pairs can 
be realized from the ground state of strongly attractive 
Hubbard model by a sudden switch of interaction from 
U < to U > 0. To sec it clearly, we calculate the 
transition probability from the ground state of attractive 
Hubbard model (marked by the symbol of cross in Fig. 
3b) to the highest state in the top band of the repul- 
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FIG. 6: (Color online) Energy vs U for states with different 
total spin S. 



sive Hubbard model (marked by the symbol of cross in 
Fig. 3a). As shown in Fig. 0b, when both the attractive 
interaction of the initial state and the repulsive interac- 
tion of the final state are strong enough, the transition 
probability is very close to 1. By this way, we can real- 
ize the repulsively paired state for the repulsive Hubbard 
model. We note that such a repulsively paired state is 
a very highly excited state with zero total momentum, 
which is different from the ^-pairing state discussed in 
Ref.QUlt. 

Next we calculate the number of doubly occupied sites 
Nd = J2i i^if^ii) f° r the FSTG state, which can be ob- 
tained by differentiation of the energy. As shown in Fig. 
[5^,, Nd decreases monotonically with increasing \U\ and 
tends to zero for U — > — oo. This indicates that there 
is no paring for the FSTG state even in the strongly at- 
tractive limit. As a comparison, we also calculate Nd for 
the ground state of the attractive Fermi gas in Fig. [5p. 
Instead, Nd for the ground state of the attractive Fermi 
gas monotonically increases to N/2 as the ground state 
is composed of N/2 pairs of fermions with the bounding 
energy proportional to U. For a bosonic STG gas, it is 
known that the STG state has even stronger local corre- 
lation than the repulsive TG gas [1, . To see whether 
FSTG state in the optical lattice has similar properties, 
we also calculate the local correlation function for the re- 
pulsive ground state. Here we note that the local corre- 
lation function (ri^n^) is noting else but Nd defined 
above. As shown in in Fig. [5£i, Nd for the ground state 
of repulsive Fermi gas also monotonically decreases to 
zero with the increase in the repulsion strength U. Nev- 
ertheless, the local correlation functions shown in Fig. 
[5^ indicate that the FSTG state has stronger local cor- 
relation than the corresponding repulsive ground state, 
which is similar to its bosonic correspondence. 

In contrast to the spinless Bose system, the ground 



FIG. 7: (Color online) Full spectra of energies vs. total 
momentum K for the HM with L = 10, N = 10, N t = TVj. = 
5, U = 15 (left) and U = -15 (right). The left inset is 
the enlargement of the lowest band of the half-filled repulsive 
Hubbard model whereas the right inset is the enlargement of 
the top band of the half-filled attractive Hubbard model. 



state of a spin-1/2 Hubbard model is highly degenerate 
in the limit of U — >• oo due to the existence of the spin 
degree j2(|. Nevertheless, the degeneracy is broken for a 
large but finite interaction strength and the true ground 
state of the repulsive fermions is the state with the low- 
est S [27]. In Fig|6l we show the ground state energy of 
repulsive Fermi gases and the lowest energy of the FSTG 
state for systems with a fixed N = 10 but various S. It 
is clear that E(S±) < E(S 2 ) for Si < ^2 on the repul- 
sive side. The energy difference between different spin 
states vanishes as U — > 00. On the other hand, for the 
FSTG state, the state with the larger S has lower en- 
ergy, i.e., E(Si) > -E(5 2 ) for Si < 52 on the attractive 
side, as shown in Fig|6l This implies that the ferromag- 
netic state with maximum S = N/2 has lowest ener- 
gies for the FSTG states with large but finite interaction 
strength. Nevertheless, energies for states with different 
total spins approach the same limit of the polarized state 
as \U\ —> co. 

Finally, we discuss the half-filling case, for which the 
ground state is a Mott state for any finite repulsion 
20]. Similar to the low-density case, we display the full 
energy-momentum spectra of the the half-filling Hubbard 
model with L = 10, N = 10, iV t = N± = 5, and U = ±15 
in Fig. [7] Similarly, the spectra is spilt into a series of 
separated bands. For the repulsive case, the lowest band 
corresponds to the Mott states with each site occupied 
by a single fermion. Above the Mott band, the sepa- 
rated middle bands correspond to the scattering contin- 
uum composed of paired fermions and unpaired fermions, 
whereas the top band is the scattering continuum of N/2 
paired fermions. In contrast to Fig. 3, the lowest Mott 
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band is obviously very narrow as the hopping process of 
a single fermion to its neighboring sites is suppressed. 
In the large U limit with U/t ^> 1, the double occupied 
states have much large energy than the states with no 
double occupancy, and the effective Hamiltonian is given 
by [H 




where Sf = ^cj a a" a >Ci, a > (a — x,y,z) are the usual 
spin operators with a a being the Pauli matrices. As 
shown in the inset of Fig. 7a, the enlarged spectrum of 
the lowest band is consistent with the spectrum of anti- 
ferromagnetic (AFM) Heisenberg model of Eq.(j8]). For 
the attractive case with U — —15, as shown in Fig. 7b, 
the spectrum has similar structure as the repulsive case 
but in a reverse order with the top band corresponding 
to the Mott states. By projecting the original Hamilto- 
nian (1) into the Hilbert spaces without any double oc- 
cupancy, we can also get the effective Hamiltonian given 
by Eq.©. Although the effective Hamiltonian has the 
same form for both the repulsive and attractive Hub- 
bard model, the effective coupling strength J has differ- 
ent sign, i.e., J = ±4£ for U > or U < 0. For U > 0, 
the effective model is an AFM Heisenberg model which 
describes the lowest continuum band of the original Hub- 
bard model. On the other hand, the effective model is a 
ferromagnetic (FM) Heisenberg model for U < 0, which 
describes the highest continuum band of the attractive 
Hubbard model. 

Given the initial state as the ground state of half-filled 
repulsive Hubbard model (or effectively the ground state 
of AFM Heisenberg model of Eq. flSJ labelled by the star 
in the inset of Fig. 7a), after a sudden switch of interac- 
tion from U = 15 to U = — 15, the state is transferred 
to the highest state of the attractive Hubbard model (or 



effectively the highest excited state of the FM Heisenberg 
model labelled by the star in the inset of Fig. 7b). By this 
way, we can effectively prepare the highest excited state 
of a FM Heisenberg model. Here we indicate the differ- 
ence from the low-density case: the final state obtained 
by sudden switch is on the top of the top band in Fig. 7b, 
whereas the FSTG state in Fig. 3b is on the bottom of 
the top band. On the other hand, if the initial state is 
the ground state of half-filled attractive Hubbard model, 
we can access the repulsively paired state in the top band 
of the repulsive Hubbard model by the sudden switch of 
the interaction to the repulsive side. 



IV. SUMMARY 

In summary, we study the properties of the FSTG 
state of the attractive Hubbard model, which is a highly 
excited state of the attractive Hubbard model without 
paired states and corresponds to the lowest real solution 
of the Bethe-ansatz equations with U < 0. Starting from 
the ground state of strongly repulsive spin- 1/2 fermion in 
ID deep optical lattices, such a state can be realized via 
a sudden switch of the interaction to the strongly attrac- 
tive regime. By calculating the transition probabilities, 
we have shown that the excited FSTG state can be ef- 
ficiently achieved in the strongly interacting regime and 
thus is possible to be realized experimentally with cold 
fermionic atoms in optical lattices. 
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